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ABSTRACT 



Although the analysis of observational data indicates 
that quasi-separatrix layers (QSLs) have to play an 
im portant role in magnetic configurations of solar 



fla fes, the corresponding theory is only at an initial 
stage so far. In particular, there is still a need in a 
proper definition of QSL. This problem is analyzed 
here on the basis of geometrical properties of the 
mapping produced by the field lines which connect 
photospheric areas of positive and negative magnetic 
polarities of active regions. In general, one can find 
on the photosphere a unique pair of locally perpen- 
dicular line elements at one footpoint of a given field 
line, so that this pair is mapped into a similar pair 
of elements at the other footpoint. Along the direc- 
tions of these elements the field line mapping only 
stretches and compresses the corresponding displace- 
ments of the neighboring field-line footpoints. This 
fact enables us to get a remarkably clear and concise 
definition of QSL as a volume filled by coronal field 
lines for which the ratio of the corresponding stretch- 
ing and compressing coefficients is anomalously large. 
The new definition is also compared with the ones 
previously introduced by other authors. 

The theory is applied to flare events of the so-called 
sigmoid-type, using an analytical model of a twis ted 
force-free configuration ( Titov fc Demoulin 1999a| , b). 
It is shown that such a configuration may contain a 
QSL even if genuine separatrix surfaces are absent. 
Identifying QSLs is an essential prerequisite for un- 
derstanding the mechanism of magnetic energy re- 
lease in this type of flares. It is also demonstrated 
that the magnetic field under study has a Hamilto- 
nian structure, which makes it possible to reveal a 
geometrical reason for the appearance of QSL in this 
case. 



ways show a relation of solar flares with the sepa- 
ratrices cau sed by magnetic null points (Demoulin 
et al. 1994| ). There is, however, a more systematic 



spatial correlation between the locations of energy 
release in flares and the regions of strong variation 
of the field line connectivity (e.g. Mandrini ct al. 
1995|, pemoulin ct al. 1997| ). Such regions, called 



quasi-separatrix layers (C^bLs) , are thought to be the 
plaus ible places for the magn etic reconnection pro- 
cess (Priest & Demoulin 1995). 



In most of the coronal volume the quasi-static condi- 
tions are fulfilled, so that the magnetic field evolves 
through a sequence of force-free equilibria. These 
conditions, however, may easily break down in QSLs, 
where due to a strong variation of the field line con- 
nectivity the rearrangement of the field lines during 
the evolution of the configuration may occur faster 
than in other places. This in turn implies a lo- 
cally large acceleration of plasma and hence a locally 
unbalanced and enhanced Lorentz force with corre- 
sponding concentrations of the current density. So 
ultimately the inertia of plasma may cause the for- 
mation of current layers in QSLs. The importance of 
inertia in the current layers at the QSLs is also fol- 
lowed from exact solutions of linearized MHD equa- 
tions describing a quasi- static evolution of inhom o- 
geneous magnetic fields (Inverarity & Titov 1997). 



The paper is organized as follows. In Sec. 2, we dis- 
cuss the difference between separatrix surfaces and 
QSLs together with the primary definition of QSLs. 
In Sec. 3, the local geometrical properties of the mag- 
netic connectivity are studied. The new geometrical 
measure for QSLs is derived and compared with the 
old ones. Sec. 4 contains an application of this theory 
to a model of flaring twisted configurations possess- 
ing a Hamiltonian structure, which plays an impor- 
tant role in the formation of QSLs. The conclusions 
are summerized in Sec. 5. 



Key words: Quasi-separatrix layers; magnetic topol- 
ogy; solar flares. 



2. THE QSL AND ITS PRIMARY DEFINITION 



1. INTRODUCTION 



Investigations of coronal magnetic fields extrapo- 
lated from photospheric magnetograms do not al- 



The magnetic field lines in solar active regions nor- 
mally connect domains of positive and negative po- 
larity of the photospheric plane, say, z — 0. Let 
the location of their footpoints in this plane be rep- 
resented depending on the polarity by the radius- 
vector r + = (x_|_,y_|_) or r = (x_,y_). The con- 
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nections of the footpoints by the field lines deter- 
mine two mutually inverse mappings LI : r + i— ► r 
and II + : r i— > r+. We shall simply use II if we 
refer to aspects valid for both mappings. Also the 
functional forms (X_(r+), Y_(r + )) = LI (r+) and 
(X_|_(r_), Y+(r_)) = JI + (r_) will be used further for 
the mappings. 

The mapping II is discontinuous at the footpoints of 
the field lines threading magnetic nulls in the corona 
or touching the photosphere, since the magnetic flux 
tubes enclosing such field lines are splitted at the 
nulls or at the touching points ( Seehafer 1986 ). The 
corresponding discontinuities serve as indicators for 
the separatrix field lines and surfaces. It is worth to 
emphasize that the coordinates (x±,y±) in this case 
need not to be Cartesian because the discontinuities 
are revealed in any system of coordinates irrespective 
of the metrics. 

However, with the help of the metrics or Cartesian co- 
ordinates one can determine not only the genuine sep- 
aratrices but also the QSLs. The integrity of the flux 
tubes is preserved within the QSLs and so the map- 
ping II remains continuous at the corresponding foot- 
points, but the shape of their cross-sections strongly 
change along the flux tubes. Thus, instead of the true 
discontinuities in II at the intersection of the genuine 
separatrices with the photosphere, there are contin- 
uous but rapid variations in II at the photospheric 
cross-sections of QSLs. These variations can only be 
detected by using the metrics, which enables us to 
measure and compare the distances between the foot- 
points in one polarity and corresponding footpoints 
in the other polarity. In this respect QSLs and sep- 
aratrices are qualitatively different objects. Indeed, 
ignoring the above metrical information about II and 
using a proper continuous change of coordinates, it 
is possible to eliminate the rapid variations in II and 
thereby the QSLs themselves, while the discontinu- 
ities of II and hence the corresponding separatrices 
are not removable in this way. However, for physi- 
cal processes with characteristic length scales much 
less than the QSL thickness such a QSL must be as 
important as a genuine separatrix. 



Fo r the determination of the QSLs Priest & Dcmoulin 
(1195) proposed to use the function (called "the 
norm ) N(r + ) or 7V(r_), which in Cartesian coor- 
dinates are 



N(t ± ) = 



dx± 
\dx± 



dy± 

dY T 
dy± 



1/2 



= 7V 4 



(1) 



It was suggested that N(r±) 3> 1 at the footpoints of 
the field lines belonging to QSLs and N(r±) « 1 oth- 
erwise. Yet this norm in application to different foot- 
points of the same field line yields generally different 
values N + and AT_, which leads to an ambiguity in 
the determination of QSLs. This disadvantage of the 
norm indicates that the adequate measure for QSLs 
must be invariant to the choice of the mapping LI 
or II, . Below we find such a measure by analyzing 
geometrical properties of the field line connectivity. 



3. REVISED DEFINITION OF QSL 



The mapping LI or LT + is locally described by its 
differential JI * or LI + * , respectively, which is a linear 
mapping from the plane tangent to the photosphere 
at one footpoint to a similar plane at the other foot- 
point. These differentials are represented by the cor- 
responding, mutually inverse, Jacobian matrices 



V 





0X_ 


dx j_ 


9y+ 






dx 4. 


9y+ 



a b 
c d 



and 



A, 



/ ax+ 

I dx- 

\ dx- 
ad — be : 



dX+_ \ 

dy- 
0Y+ 

dy- ) 
= det p_ . 



d -b 
—c a 



(2) 

. (3) 
(4) 



We assume hereafter that (x± , y± ) are measured in 
one Cartesian system of coordinates covering the 
whole photospheric plane. Eqs. (g) and (^) show 
that it is sufficient to have only one of these matrices 
for a local description of the magnetic connectivity. 

An orthonormal vector basis (u+,v + ) in the tan- 
gent plane at some footpoint r + is generally mapped 
by II * into a non-orthonormal basis (u_ , v ) with 
u = J2u+ and v = + 22v + in the correspond- 
ing tangent plane at Jl (r+) (Figure [l]a). One can 
choose, however, a special orientation of (u + ,v + ) 
such that the basis (u_ , v_ ) becomes orthogonal 
(Figure jl|b) . Indeed, let the angles 7+ and 7_ , re- 
spectively, determine the directions of u+ and u_ 
with respect to the x-axis of the global Cartesian 
system of coordinates. 7?. 7+ and 7£ 7 _ are the cor- 



responding matrices of rotations, 
from Figure |l]b that |u_| T?.^ 1 u+ : 



It follows then 
and 



n~l p_ u+ 



n- ; v 



' " — TZ-^ v+ or, in other words, 



diag(|u_|,|v_|)=7£- 



(5) 



It means that the two non-diagonal elements at the 
right hand side of (||) must vanish. This condition 
can be satisfied by 



7+ - 7- 
7+ +7- 



= arg£, 
= arg C, 



(6) 
(7) 



where the function arg determines the arguments of 
the complex values (i = y— T) 



i = a + d + i(b-c) } 
C = a-d + i(b + c). 



(8) 
(9) 



This solution for 7+ and 7_ can be checked by sub- 
stitution of (|^) and the corresponding expressions for 
the rotation matrices in (0) . Other solutions are also 
possible but not of interest, since they are simply the 
result of mirror reflections about the directions given 
by (| and |). 



Eqs. 



and (@) yield 



|u_| = (|ei + ICI)A 

iv_i = (lei-ici)A 



(10) 

(11) 




Figure 1. The mapping of an orthonormal basis (u+,v+) by the differential II* (a,b): its image basis (u_,v_) is in 
general non- orthonormal (a), but it is orthogonal for a special orientation of (u+, v+) (b) when |u_ + v_ | determines 
the norm iV(r+). The inverse differential II* maps the orthonormal basis (u_,v_), rotated with respect to (u_,v_) 
on 7r/2, into the orthogonal basis (u+, v+) (c), so that |u+|/|v+| = |u_ |/|v_ | but |u+ + v + 1 = N( II (r+)) 7^ N(r+). 
The shaded areas show how these properties of the mappings manifest in the corresponding cross-sections of thin flux 
tubes. 



which means that |u_| > |v_| in the solution under 
consideration. By using these expressions and 
(^) one can obtain 

|u_ + v_| = (a 2 + b 2 + c 2 + d 2 ) 1 ' 2 = N+. (12) 

Thus the norm N + determines simply the length of 
the rectangular diagonal constructed on the vectors 
u_ and v_ (Figure ^jb). 

It is worth to determine similar characteristics for 



the reverse differential JI + * at r_ = + n (r + ) . The 
simplest way to do this is just to change in (p|)"(|ll) 
the superscripts ± on =p and the elements of^D (|5 ) 
to the corresponding elements of P + (||). Then we 
obtain for the complex values similar to (||) and (Q) 

£ = [d + a-i(b-c)]/A + = ^/A + , (13) 
C = [d-o-i(6 + c)]/A + = -C/A +> (14) 
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yielding the following equations for the angles: 

7- -7+ = arg£=-arg£, (15) 
7- +7+ = argC + 7r. (16) 

They show that the orthonormal basis (u_, v_) cor- 
responding to this solution is turned on 7r/2 with re- 
spect to the basis (u_ , v ) and the same is valid for 
(u+. v+) and (u+,v+) (Figure Jljb, c). In application 
to jlQ ) and (|ll]) the consideredtransformation gives 



Iv 4 



= l/|v_| 
= l/|u_| 



(17) 
(18) 



So the norm in the negative polarity at the footpoint 
r_l= JI (t + ) is 



-ftz 



JL = JV+/ | A+ | 



This consideration shows that the determination of 
QSLs by means of the norms N + and A_ must really 
lead to different results if |A + | ^ 1. Combining of 
these results in order to find the "proper" QSL is in 
fact equivalent to finding only in the positive polar- 
ity the regions where the functions N+ and N+/\A+\ 
acquire anomalously large values. Yet the simultane- 
ous using of two different functions for characterizing 
QSLs does not look as a theoretically well-founded 
approach. 

One can resolve this difficulty by noticing, first, that 
the mapping II can locally be described by |u_|, |v_|, 
7+ and 7_, where only |u_| and |v_| determine the 
value of footpoint displacements, while 7+ and 7_ 
define their directions. So it would be natural if the 
required characteristics is a function of |u_| and |v_| 
only. We think that the ratio |u_ |/|v_ | > 1 is a more 
suitable quantity for describing QSLs than the con- 
sidered above norm. According to ( |i~7j ) and ( |l8| ) (see 
also Figure |l]b, c) it coincides with the ratio |u~[~[/|v + | 
and thereby characterizes the magnetic connectivity 
itself rather than one of the mappings JI or JI + . It 
determines the degree of flattening of the elementary 
flux tubes at their photospheric ends (see shaded re- 
gions in Figure P),c), so the tubes characterized by 
extremely high values of this Quantity should be re- 
ferred to as QSL. By using (p|)— (|TT|) one can derive 
that 



u_ 



|v_| = Q/2 + VQ74-1, (20) 
Q = Nl/\A + \, (21) 



which shows that |uj/|v_| » Q for Q 3> 1. It is 
seen from ( |l9| ) and (2W that Q is simply a product 
of the norms N + and~7V_ calculated at the opposite 
ends of the field lines. Geometrically, the value Q is 
the ratio of the square of quadrat with the side N + 
to the square of the rectangular based on the vec- 



tors u and v_, because Aj 



|v I. Note 



also that Q > 2, since N + is the length of the diag- 
onal in this rectangular. The expression for Q has 
an elegant form, invariant to the change + <-> — 
and simpler than (EOI), so it is reasonable to define: 



The QSL is a flux tube consisting of magnetic field lines 
with Q > 1. 

One can see from (|2l]) that if the determinant A + 
varies in the photospheric plane as strong as N + , the 
QSLs defined with the help of Q or N + have to be 



different, otherwise they must be approximately the 
same. 

Some comments should also be added on another 
characteristics for determining QSLs, namely, on the 
so-called differential flux volume 

f+ n (r+) d/ 
V(r + )= I -^V(Jl(r + )) 1 

in which the integration is carried out along the cor- 
responding field line. This value has appeared in the 
analysis of current sheet formations along separatrix 
surfaces in quasi-static evolutions of 2 gD magnetic 



configura ti ons ( Zwingmann et al. 1985 , Low fc Wolf- 



son 1 988. Vekstein et al. 1991, Vekstein & Priest 



1992 ). It has also been used (under the name "delay 



function" ) for studying 3D mag netic topo logy caused 
by the presence of null p oints ( Lau 1995 ). Recently 
Bchindler fc Birn (1999) have shown that strong spa- 
tial variations of V may cause a similar formation 
of current layers in 3D magnetic fields. So it seems 
natural at first sight to use this value for character- 
izing QSLs. However, if one fixes the footpoints of 
the magnetic field lines in a given configuration and 
exposes the coronal volume to a smooth deformation, 
then V will be changed depending on such a defor- 
mation, while the field line connectivity will remain 
the same. Thus, the quantity Q as a measure of the 
field line connectivity is better than V, although this 
does not reduce the value of V in understanding of 
the above mentioned current sheet formation. 



4. HAMILTONIAN STRUCTURE AND QSL IN 
FLARING TWISTED CONFIGURATIONS 



The topology of a twisted magnetic configuration 
provides im portant clues for understan ding sigmoid- 
like flares ( fTitov fc Dcmoulin 19994 b )- Such a 
configuration can be modeled by a force-free circu- 
lar flux tube with total current I, which is embedded 
into potential field produced by a pair of magnetic 
charges q, —q and a line current Jo (Figure |2j). The 
charges are on a distance 2L from each other and lie 
together with the current Jo on the axis of symme- 
try of the configuration, which is on the depth d be- 
low the photospheric plane z = 0. Below this plane 
the sources and their fields have no a real physical 
meaning: they are used only to construct the proper 
magnetic field in the corona. The photospheric ver- 
tical component of this field and current density of 
the flux tube can be considered a posteriori as the 
corresponding boundary conditions for the resulting 
configuration. 

The construction of this magnetic field excludes the 
existence of any null point in our model, so the topol- 
ogy of the coronal field may be non-trivial only due to 
the presence of t he so-called bald patches (BPs, see 
Litov et al. 1993), which are the segments of the pho- 
tosplicric inversion line (IL) where the coro nal field 
lines touch the photosphere ( |Seehafer 1986 ). Thus, 
the separatrices, if present, consist of only the field 
lines starting at the BPs. The full topological analy- 
sis of the field then reduces to determining the BPs 
and their associat ed separatrices. Such an a nalysis 
has been made by Litov & Dcmoulin (1999) for the 
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scenario in which the major radius R of the flux tube 
(Figure grows, starting from R = d, while the 
intensity and positions of sub-photosphcric sources 
remain constant. It was found that the BP appears 
first at some R = i? a > d and then grows in size with 
increasing R. After reaching some R — R\> ^> R a 
the BP bifurcates at the center of the configuration 
to give birth to a generalized separator field line at 
the intersection of two separatrix surfaces. The gap 
between bifurcated parts of BP grows with R, while 
these parts shrink and eventually disappear at some 
R = i?d (> Rb), where a transition to a topologically 
simple arcade-like configuration occurs. 




Figure 2. The scheme of the model for the twisted mag- 
netic configuration 

In spite of this simplicity, the resulting magnetic field 
might have QSLs, which are as important as the gen- 
uine separatrices under the appropriate conditions. 
The aim of the further consideration is to check this 
possibility by using the above theory of QSLs for the 
same scenario of the magnetic field evolution. 

Here we use the cylindrical system of coordinates 
(r±,6,x), whose axis is directed along Iq and where 
r± = \r±\ with r_L = (0,y, z + d), 9 is the angle 

between rj_ and y, and 9 is the corresponding unit 
vector. Due to the axial symmetry about the line cur- 
rent Iq (Figure ||) the magnetic field can be written 

as 



B 



v x [ — e 

. r± 



Bn 9 . 



which, in particular, gives 

1 a* 



B x — 



9rj 



(22) 



(23) 



It is convenient in our case to use 9 for parameterizing 
field lines, so that their equation dr x B = with 
dr = dr± r± + r± d9 9 + dx x yield 



da; 

d9 
dr^ 
~d9~ 



1 



a* 

i drj_ ' 
1 9* 
Be ~dx ' 



(24) 
(25) 



These are the components of the field line equation 
along fj_ and x, respectively, while the component 



along 9 requires d^ = al ong 
is simply a consequence of ( p4|) 



the field lines, which 

-(H). 



The axial symmetry allows also to write 
component of the current density j as 



the 



1 



Jx = 



d 



H r± drj_ 



from where after introducing the whole current J = 
f Q ± j x 2irr '_l dr± through the perpendicular to the x- 
axis disc of radius r±_ one can obtain 



Bn 



/WW 

2nr I 



(26) 



Here we put J = J(^), since the current lines lie on 
magnetic surfaces ^ = const (note that in our model 
j || B and B • V* = according to Eq. More- 
over, this is also valid in more general case of mag- 
netostatic configurations. By using (26) and normal- 
izing the coordinates as f± = rj_/L, x — x/L with 
subsequent omitting of the tilde, one can transform 
(pi|)-(p5|) to the following one-dimensional and au- 
tonomous Hamiltonian system: 



dx 
dff 

dm 

~dJ 



dH 
dvj ' 
OH 

dx 



(27) 
(28) 



where 9 is a time-like variable and instead of rj_ the 
new momentum-like variable 

w = ]nr± (29) 

is introduced together with the Hamiltonian function 



H 



2n 



d* 



J J(*) 



(30) 



The Hamiltonian structure is revealed here only with 
the help of the axial symmetry and the magnetostatic 
condition J = J^), so it must be inherent to both 
force-free and magnetostatic fields having the axial 
symmetr y. The analog ous fact has been early estab- 
lished by Lewis (1990) for toroidal geometry but we 
derived it m a much simpler way by using a different 
Hamiltonian and canonical variables. 

This structure provides a powerful tool for investi- 
gating the geometry of magnetic field lines in our 
model. Each field line can be considered now as a 
trajectory in the extended phase space {x, vj, 9) of 
the autonomous one-dimensional dynamical system 
(p7|), (p8|). In the following it turns out that the 
most interesting behavior of the field lines is outside 
of the flux tube, where the magnetic field is purely 
potential and J(^) = Iq, so that 



Bo = 



The Hamiltonian 



H 



2frr± 
Of) then reduces to 
2tt<I' 



(31) 



(32) 
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which is simply the poloidal flux function multiplied 
by a constant. This function is the following linear 
superposition: 



(33) 



where the flux functions ^± q and respectively, 
correspond to the charges ±q and flux tube. The 
x-component of the field produced by one charge is 



B±„ x = 



±q(x T 1) 



r 2 + ( x =p l)2]3/2 ' 



which together with (E3L) enables us to obtain 



(34) 



(35) 



Let the minor radius a of the flux tube be much less 
than its major radius R, then ^ j must be well ap- 
proximated outside the tube by the flux function of 
the circle coil of radius R with current /. Due to 



the axial symmetry tyj(r±,x) = rj_ Aj (r± , x) , where 
A{ is one non-vanishing y-component ot the vector- 
potential 



Aj(r) = 



Ho I 

47T 



dr' 



r - r' 



in which the integration is made over the whole coil. 
After some transformations it gives 



Ho I 

V8 



where 



{LRr±v) 1/2 ^(i 



2R/Lr± 
R 2 /L 2 + r 2 ± +x 2 



$(w) 



cos(9d<9 



(1-VCOB9) 1 / 2 



One can verify that $(u) is approximated by 
$(u) « 0.043 u - 0.089 u 2 - 0.207 ln(l - v) 



(36) 

(37) 
(38) 

(39) 



with the relative error < 1% on the interval [0,0.9], 
which is good enough for our further purposes. 

The Lorenz force F q = -2qLI(R 2 + L 2 )- 3 / 2 caused 
by the interaction of the current / with the fields 
B± qx and the self-force Fj — ^o/ 2 (ln8i?/a — 
5/4)/ (47ri?) of the flux tube due to its curvature coun- 
terbalance one another if the current is 



8TrqLR{R 2 + L 2 )- 3 ' 2 
fi [H8R/a) - 5/4] 



(40) 



The expression for Fj used here corresponds to the 
simplest case when the equilibrium current I in the 
tube is uniformly distributed over its circular cross- 
section of radius a. Other possible current dis tribu- 
tio ns iust slightly modify the constant 5 /4 (see [Titov| 
& fjemoulin 1999] ). For the assumed small radius of 
the flux tube (a <c R) the toroidal field Bg is ap- 
proximately homogeneous in the vicinity of the tube. 



This enables us to relate the radius a to the total 
number nt of field line turns in the tube as follows: 



HqIR 
2tt a 2 \B 



IR 2 
JloV 2 



(41) 



In the above mentioned scenario of the emerging flux 
tube only the equilibrium current / and minor radius 
a change with growing major radius R according to 
(Eo|) and (pd|), while n t and the other parameters of 
the model q, d, L, and Iq remain constant. At each 
moment of such an evolution Eqs. (p^), (|3^), (|3^), 
and ( p3q) -( |4l|) describe the Hamiltonian dynamic sys- 
tem (27])- ipq) corresponding to the magnetic field 
outside the nux tube. 

The Hamiltonian in this region has a hyperbolic 
critical point, which is in fact the X-point of the 
poloidal magnetic field. Figure y shows the corre- 
sponding phase portrait of the system for the follow- 
ing set of parameters: nt = 5, q = 100 T • Mm 2 



L = 50 Mm. R = 110 Mm, I = -7 TA (cf. |Titov 
fc Dcmoulin 1999| b) . This portrait is a result of pro- 
jection (x,w,0) — > (x,xu) of several magnetic field 
lines. It can also be considered as a set of contours 
H = const, since our autonomous phase flow con- 
serves the Hamiltonian. Taking into account Eq. ( [32] ) 
or more general Eq. ( |30j ) , one can obtain another use- 
ful interpretation of the portrait, namely, as a cross- 
section 9 = const of the corresponding poloidal mag- 
netic surfaces \1/ = const shown in the coordinates 
(x, w). These coordinates differ from Cartesian ones 
only by the logarithmic transformation (w = lnrj_) 
of the vertical axis, so Figure || yields a rather real- 
istic information about the shape of these surfaces. 



1 

0.8: 

0.6 : 
03 0.4 
0.2 


-0.2 
-0.4 




-0.6 -0.4 -0.2 0.2 0.4 0.6 
X 



Figure 3. The phase volume at four different equidistant 
moments: the initial volume is represented by a small cir- 
cular, which can be considered as the corresponding pho- 
tospheric cross-section of the flux tube starting near the 
center of the configuration; the final highly stretched vol- 
ume corresponds to the moment when the field line, pass- 
ing through the X-point of the poloidal field, reaches the 
photosphere. 
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The important property of autonomous Hamilto- 
nian systems is the conservation of the correspond- 
ing phase volume in the phase flow. In our case this 
volume is two-dimensional, so the area of the plane 
(x, vj) is conserved in such a flow. This property in 
combination with the presence of the hyperbolic crit- 
ical point, provides a mechanism of the QSL forma- 
tion in our configuration. Indeed, any cross-section 
of an elementary flux tube corresponds to an area in 
the phase plane and the variation of this cross-section 
along such a tube corresponds to an evolution of the 
appropriate phase area due to the phase flow. Given 
a small circle near one of the separatrices on the phase 
portrait as an initial area, its area-conserving evolu- 
tion with the divergent character of the phase flow 
near the critical point will inevitably cause a stretch- 
ing and compressing of this area across and along, 
respectively, the direction of the flow (see Figure 0). 
This means that the elementary flux tubes starting 
nearby the poloidal separatrix surfaces have to expe- 
rience especially strong flattening, which is a charac- 
teristic feature of the field lines in QSL. 

One can also expect from the above consideration 
that the maximal value of Q = Qx is achieved at 
the field line passing through the X-point of poloidal 
magnetic field. It is possible to derive an analytical 
expression for Qx by first linearizing Eqs. (p7|)-(p8|) 
at the hyperbolic critical point of H and then solving 
them in order to calculate the corresponding deriva- 
tives for Qx- The final result reads 

Qx = 2+ f4+ f 9 V rx , 2 2 ) sinh 2 (AA#) , (42) 
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Figure 4- The dependence of Q = Qx ° n the radius of 
the flux tube R for the field line passing through the X- 
point of poloidal magnetic field at different values of the 
line current Iq (Figure ^j. 
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where rx ( < R/L) stands for the dimcnsionless ra- 
dius r± corresponding to the X-point of the poloidal 
field. The value A is given by 



A = rx 



d 2 H 



dx 2 



(43) 



x— 0,t37— In rx 



and according to (|22|), (|3~T| ) and (32) it characterizes 
the strength of the pofolaal fieldon the distances 
~ Lrx in comparison with the toroidal field at the 
X-point. The value 



A6 = ir — 2 arcsin 



Lr x 



(44) 



determines the arc length of the coronal part of the 
field line passing through the X-point. 

Eq. ( |42] ) shows that Qx grows exponentially with 
A AO, so Qx strongly depends on the relationship be- 
tween poloidal and toroidal components in the vicin- 
ity of the X-point and the length of the field line pass- 
ing through this point. This value has been computed 
for the above mentioned parameters and d = 50 Mm 
in the scenario of the emerging flux tube. The de- 
pendence Qx(R) at several different values of Jo (see 
Figure (^)) suggests that the QSL has to appear in 
the twisted configuration when the major radius of 
the flux tube reaches R <; 100 Mm. This corresponds 
approximately to the moment when the gen uine sep- 



aratrix surfaces (caused by BPs) disappear (Titov 
Debioulin 1999|b). 



In this paper, new results concerning the theory of 
the quasi-separatrix layers (QSLs) in coronal mag- 
netic fields are reported . We have revised the pre vi- 
ous definition of QSLs ( Priest fc Demoulin 1995| ) by 
analyzing the local geometrical properties of the field 
line connectivity. This analysis provides a new geo- 
metrical measure Q, which determines the degree of 
flattening of the elementary flux tubes. Co ntrary to 
the so-called norm N proposed previously by Priest & 
Demoulin (1995), the new measure Q yields the same 
results tor both positive and negative magnetic po- 
larities on the photosphere. In other words, Q char- 
acterizes the magnetic connectivity itself rather than 
the corresponding field line mappings (from positive 
to negative polarity or vice versa). 

The old and new measures are related by a very sim- 
ple formula Q = N 2 / A, in which A is the Jacobian 
of the corresponding field line mapping. Thus, if A 
varies on the photosphere as strong as N, the new and 
old measure must yield different locations for QSLs. 

We have applied this theory to the model of 
twisted magnetic configuration proposed in ([Titov 



Demoulin 1999ap) for the analysis of magnetic topol- 



ogy in sigmoid-likc solar flares. It has been demon- 
strated that the appearance of QSLs in such a model 
is related with the Hamiltonian structure of the mod- 
eling magnetic field: the conservation of the phase 
area in the corresponding phase flow and the pres- 
ence of the X-point in the poloidal magnetic field 
provides the mechanism for the formation of QSLs 
in this configuration. The results of computations 
of Q suggest that the QSL must be formed in such 
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a configuration approximately after that the genuine 
separatrix surfaces disappear when the twisted flux 
tube emerges high enough into the corona. 
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